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DESCRIPTION 

METHOD AND DEVICE FOR THERMAL THROUGHFLOW MEASUREMENT 
WITH PULSED HEATING POWER 

TECHNICAL FIELD 

The present invention relates to the field of measuring 
flows of gas or fluid using thermal sensors. It is 
based on a method and a sensor for mass flow measure- 
ment according to the preamble of the independent 
claims . 



PRIOR ART 



US Patent No. 4,501,145 discloses a generic method and 
device for thermal flow rate or flow velocity 
measurement. Here a sensor is heated with a heating 
pulse of constant power, at least one build-up time is 
measured until a preset temperature threshold is 
reached and from this a fluid parameter, such as e.g. 
the flow rate, dependent on the heat transfer 
coefficient is measured. The disadvantage is that the 
rise in temperature proves to be more and more 
undifferentiated with increasing flow rate and comes to 
saturation, so that adequate measuring sensitivity is 
achieved only in a restricted measuring range of flow 
rates. In addition, the heating power or heating power 
and the temperature threshold are to be carefully 
matched to one another and to the desired flow rate 
measuring range . 

EP 0 180 974 discloses a process and a device for 
measuring flow rates or mass flows. Thereby, flow- 
dependent sets of characteristic curves are determined 
for the unstationary probe temperature evolution at a 
constant heating power and are related to the flow 



- 2 - 



rate, in that temperatures are measured at preset times 
or time intervals until predetermined temperature 
threshold values are reached. The disadvantage is that 
the flow rate measuring range is restricted to the 
evaluable ranges of the sets of characteristic curves 
on account of the again decreasing temperature 
resolution at high flow rates. 

DESCRIPTION OF THE INVENTION 

The object of the present invention is to provide a 
method and a device for pulsed mass flow measurement 
with an improved measuring sensitivity in an enlarged 
measuring range. This object is achieved according to 
the invention by the features of the independent 
claims . 

In a first aspect the invention relates to a method for 
measuring a flow rate or a mass flow of a fluid, in 
particular for measuring hot water supply in the 
private, public or industrial sector, in which the 
fluid is guided over a sensor element, which has a 
heating means for inducing temperature changes and a 
sensor means for determining its temperature, wherein 
at least from time to time the heating means is 
operated with a heating power in the form of heating 
pulses and a flow-dependent threshold value time is 
measured at the sensor means until a preset temperature 
threshold value is reached, wherein during at least 
some of the heating pulses a non-constant heating power 
with a substantially sublinear build-up dynamics as a 
function of time is selected in order to at least 
partially compensate a nonlinear behaviour of the 
threshold value time as a function of the flow rate or 
flow velocity. A sublinear build-up dynamics means that 
during a heating pul se the heating power is a function 
of time with a monotonously decreasing first time 
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differentiation. Substantially sublinear means that 
short-term deviations therefrom, e.g. a short-term 
increasing time differentiation, is tolerable, as long 
as the build-up dynamics remains globally, i.e. under 
formation of section-wise average values over the 
entire heating pulse, sublinear and has a first time 
differentiation decreasing from section to section. 
Such a flattening increase in heating pulse output 
leads to the desired effect, that both low and high 
flow rates can be determined with great precision and 
at the same time relatively short measuring duration. 
With the method the measuring range can therefore be 
expanded, the measuring precision can be increased, the 
measuring time can be shortened and, if required, the 
measuring repetition rates can be increased. In 
addition, it is not necessary or uncritical to adapt 
the heating power to the temperature threshold or the 
flow rates-measuring range, because the temperature 
threshold is reached in any case owing to the increase 
in heating power. 

In a first embodiment the build-up dynamics is varied 
as a function of time and optionally of the flow rate 
itself to be measured, such that the threshold value 
time is a linear function of the flow rate at least by 
approximation, i.e. in particular at least on discrete 
flow rate values. Due to the linearity a measuring 
sensitivity substantially uniform over the entire flow 
rate measuring range is reached. At the same time 
smaller local deviations from the linearity are 
tolerable. The linear measuring characteristics can be 
achieved, if required, at least pointwise over the 
entire measuring range, in that adapted build-up 
behaviours of the heating pulse power are used for 
different flow rate partial measuring ranges. 
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In the embodiment according to Claim 4 such an adapted 
build-up dynamics is achieved by the flow rate- 
dependent amplitude factor. Thereby, a linear 
characteristics of the threshold value time as a 
function of the flow rate is realised at least point- 
wise and independently of the first thermal transfer 
resistance between heating means and sensor element 
surface . 

The advantage of the embodiment according to Claim 3 is 
a particularly simple and easily implementable time 
dependence of the build-up dynamics of the heating 
power. This time dependence is suited in particular to 
the improved linearisation of the threshold value time 
as a function of the flow rate, the less the first 
thermal transfer resistance is. 

The advantage of the embodiment according to Claim 5 is 
that for a simple cylindrical form of the sensor 
element an exact flow dependency of the second thermal 
transfer resistance can be given. 

The advantage of the embodiments according to Claims 6 
and 7 is that the implicit problem of the dependency of 
the heating power on the variable to be measured is 
resolved simply and reliably by in-advance determined 
calibration curves and by estimation or a-priori 
knowledge of the presumable flow rate. 

The advantage of the embodiment according to Claim 8 is 
that with negligible first thermal transfer resistance 
the linear relation between threshold value time and 
flow rate can be calculated over the entire flow rate 
measuring range validly and exactly. 

In a second aspect the invention relates to a mass flow 
sensor for detecting a flow rate or a mass flow of a 
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fluid according to the above-described process. The 
sensor comprises a sensor element with a heating means 
and a sensor means for thermal measuring in a fluid and 
a control and evaluating processor unit with a heating 
control for generating heating pulses for the heating 
means and a measuring device for evaluating the thermal 
measurement and for determining a flow rate or a mass 
flow from a flow -dependent threshold value time until a 
preset temperature threshold value at the sensor means 
is reached, wherein the heating control comprises means 
for generating a non- constant heating power with a 
substantially sublinear build-up dynamics as a function 
of the time, and the control and evaluating processor 
unit has means for at least partial compensation of a 
nonlinear behaviour of the threshold value time as a 
function of the flow rate. 

The advantage of the embodiments according to Claims 10 
and 11 is a particularly simple and precise sensor 
control and detection of measuring variables. 

Further embodiments, advantages and applications of the 
invention will emerge from the dependent claims as well 
as from the following description and the figures. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates in cross-section a pipe flowed 
through or a flow channel with an inventive 
mass flow sensor for operation with non- 
rectangular heating pulses; 

Figures 2a, 2b, 2c, 2d illustrate pulsed operation of a 
mass flow sensor with rectangular heating 
pulses according to the prior art; 

Figures 3a, 3b, 3c illustrate pulsed operation of a 
mass flow sensor with inventive, non- 
rectangular heating pulses; and 
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Figure 4 illustrates a set of calibration curves of 

the threshold value time as a function of 
the flow rate for three flow-dependent 
heating powers. 

In the figures identical parts are designated with the 
same reference numerals. 

WAYS OF IMPLEMENTING THE INVENTION 

Figure 1 illustrates a thermal mass flow sensor 1, 2, 
which comprises a sensor element 1 arranged in a flow 
channel or pipe 4 and a control and evaluating 
processor unit 2. A fluid 3, e.g. a liquid 3 or a gas 
3, flows in the pipe 4 with a flow or rate profile 5. 
The sensor element 1 is subjected to a flow rate v to 
be measured. The sensor element 1 has a heating means 
la for inducing temperature changes and a sensor means 
lb for determining its temperature. 

According to Figures 2a to 2d, it is known to operate 
the heating means la with a heat output or heating 
power P in the form of constant heating pulses 6, to 
measure a flow- dependent temperature change T at or on 
the sensor element 1, and from this - with a given con- 
stant, i.e. rectangular heating power - to draw 
conclusions about the flow rate v or the mass flow. 
Instead of temperatures T at a fixed measuring time t, 
also flow- dependent threshold value times t s can be 
evaluated until a preset temperature threshold value T s 
is reached. Figure 2b shows by way of example two time- 
dependent temperature behaviours T(t) for two flow 
rates v a <v b . Figure 2c shows the resulting measuring 
curve or characteristics of the temperature at a 
specific time as a function of the flow rate T(v), 
which has a strongly decreasing measuring sensitivity 
for increasing flow rates v. The measuring precision 
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and the measuring range are thereby restricted in a 
very unfavourable manner. Figure 2d shows the 
corresponding strongly nonlinear behaviour of the 
thermal measurement of the threshold value time t s as a 
function of the flow rate v. The nonlinear t s (v) 
characteristics has a number of drawbacks. In the event 
of low flow rates v, the measuring sensitivity is 
minimal. At high flow rates v the theoretically high 
measuring sensitivity cannot be capitalised on, because 
the threshold value times t s can be determined only 
with major uncertainty from the flat, glancing 
intersections between the temperature rise T(t) and the 
threshold value temperature T s . Moreover, setting the 
threshold value temperature T s is highly critical. If 
T s is chosen to be low, the threshold value times t s 
become short and accordingly t s (v) becomes very flat 
and the measuring resolution for small v becomes poor. 
If T s is chosen to be high, the nonlinear, over - 
proportional increase of t s (v) already starts at low 
rates or velocities v and shows a steep incline because 
of the flat intersections. In the worst case scenario 
T s is not reached and the flow rate v can no longer be 
measured. Also in the permitted measuring range for v a 
subtle choice and, if required, repeated adapting of T s 
is required to obtain a useful measuring range with an 
evaluable measuring characteristics t s (v) . 

In figures 3a to 3c the method according to the present 
invention is explained. At the time t Q a time -dependent 
heating power signal P(t), possibly dependent on the 
selected flow rate range, is started and is increased 
according to non- rectangular , sublinear build-up 
dynamics P(t). Such heating pulses 7 are supplied to 
the heating means la at least from time to time (not 
illustrated) . Preferably, the flow-dependent times ti, 
t 2 are measured until a preset temperature threshold 
value T s is reached and the flow- dependent threshold 
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value times t s -ti-t 0 for v a and t s =t 2 -t 0 for v b are 
determined therefrom. In contrast to Figure 2c, t s (v) 
is a monotonously increasing function. On account of 
the inventive flattening of the build-up dynamics P(t) 
the interfering saturation behaviour in T(t) (Figure 
2c) or the interfering nonlinear behaviour in t s (v) 
(Figure 2d) can now extensively be compensated, 
resulting in an easily evaluable measuring function 
t s (v) that is useful over a wide measuring range of 
flow rates v. Due to the expanded measuring range the 
mass flow sensor or flow rate sensor is particularly 
suited to measuring hot water supply in the private, 
public or industrial sector. 

A detailed analysis for an optimal configuration of the 
heating power pulses 7 is given hereinbelow. In a 
simple thermal model the sensor element 1 has a heat 
capacity C s , a first thermal transfer resistance R s 
between the heating means la and a surface 10 of the 
sensor element 1 and a second thermal transfer 
resistance R I= (h«A) -1 between the surface 10 of the 
sensor element 1 and the fluid 3, wherein h is a flow- 
dependent heat transfer coefficient between the sensor v.. 
element 1 and the fluid 3 and A is a contact surface . 
between the sensor element 1 and the fluid 3. According 
to VDI heat atlas (VDI Guidelines 3522, VDI Measuring 
Handbook I, VDI-Verlag GmbH, Dusseldorf 1987), the heat 
transfer coefficient h is given, for a cylindrical, 
laterally flow-exposed sensor element 1 with diameter 
d, by the equation 

h = (X/d) • 1.11 • C • Pr 0 - 31 • Re m , (Gl) 

wherein X is a heat conductivity of the fluid 3, C is a 
parameter and m an exponent, that both depend on a 
Reynolds number Re of the fluid 3, and Pr is a Prandtl 
number of the fluid 3. With Pr=rj*c P /X, with r| being a 
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dynamic viscosity of the fluid 3 and c P being a 
specific heat of the fluid 3, and Re=p*d»v/r|, with p 
being a density and v being a flow rate or velocity of 
the fluid 3, h is proportional to v m and the following 
is valid 

Ri = Y • v~ m (G2) 

wherein y=d/ (A # A # 1. 11 *C • Pr 0 * 31 • (p •d/rj ) m ) is a constant. 
For the Reynolds number- dependent parameters C and m 
tabulated values can be used, e.g. C=0.615 and m=0.466 
for Reynolds numbers between 40 and 4000, which 
characterise a laminar flow of the fluid 3. 

In the following a numerical example is given for water 
at room temperature: rj = 1.10~ 3 Ns/m 2 , p=998.2 kg/m 3 , 
A=0.598 W/(m*K) and C P =4182 J/ (kg*K) . Then Pr=6 . 993 and 
for d=1.9 mm Re=1897*v. For rates v between 0.02 m/s 
and 2 m/s the parameters C and m assume the values 
C=0.615 and m=0.466. This leads to heat transfer 
coefficients h between 2140 W/ (m 2 *K) at v=0.02 m/s and 
18270 W/(m 2 *K) at v=2 m/s and corresponding second 
thermal transfer resistances Rx=6 . 5 K/W at v=0.02 m/s 
and R^O.7 K/W at v=2 mfs for a sensor element 1 having 
a length of 1=12 mm. 

With a time -dependent and speed- dependent build-up 
dynamics P(t) of the heating power according to the 
equation 

P(t) = P 0 • (1+Rs/Ri)" 1 • t a (G3) 

wherein P 0 is a heating power factor and a is an 
exponent of the build-up dynamics P(t), the above- 
mentioned nonlinear, in particular over-proportional 
behaviour of the measuring characteristics t s (v) can 
approximately be compensated or even the measuring 
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characteristics t s (v) can be linearised at least 
partially or pointwise at sampling points. In 
particular t s >x is valid for measuring times or 
threshold value times, in particular t s >10«T, with 
T=C S *R S 

T(t)= P 0 • t a • R x + T F , (G4) 

wherein T F designates the undisturbed fluid 
temperature. The same heating behaviour T(t) is 
obtained in the approximation R s <Ri, in particular 
R s /Ri<0.1 and particularly preferred R S /R I <0.01 / with a 
speed- independent heating pulse output P (t) = P 0 • t a . 

Assuming a choice of the heating power exponent a = m, 

T (t) = P 0 • y • (t/v) m + T F (G5) 

is valid, and the threshold value time t s for a given 

temperature threshold value T s is expressed by the 
equation 

t s (v) = (T s - T F ) 1/m • (P 0 • y)- 1/m * v . (G6) 

The measuring duration or threshold value time t s is 
also a linear function of the flow rate v to be 
determined, as shown in Figure 3c, if the heating power 
exponent a is identical to the rate exponent m of the 
heat transfer coefficient h, and, in particular, if the 
first thermal transfer resistance R s of the sensor 
element 1 is negligibly small compared to the second 
thermal transfer resistance R x to the fluid 3. The 
latter condition is satisfied if Rs/R^l, preferably 
Rs/Ri<0.1 and particularly preferred R s /Ri<0.01. 

The build-up dynmics P(t) is also selected 
advantageously proportionally to t m , wherein the 
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exponent m is chosen smaller than 1, in particular 
m<0.5 and particularly preferred m=0.466 for a Reynolds 
number between 40 and 4000. The heating power exponent 
a can also be selected identical to m only by 
approximation, i.e. oc^m, in order to achieve the 
desired, at least approximate compensation of the 
nonlinear behaviour of t s (v) at high flow rates v. 

The build-up dynamics P(t) is furthermore selected 
advantageously proportional to a time- independent 
amplitude factor (l+Rg/R!)" 1 . 

In the general sense the heating power according to 
equation (G3) is varied via Rj also in dependence on 
the flow rate to be measured or the presumed measuring 
range of the flow rate v to be measured, in such a way 
that there is again a linear relationship between the 
t threshold value time t s and the flow rate v. This 

p linear relationship is also independent from R S/ if R s 

B assumes a non-negligible value relative to R x . 

In this case the measuring procedure is preceded by a 
calibration and a selection of a matching calibration 
curve 8 according to Figure 4. In a first procedural r i 
step discrete values of the flow rate v ± are selected 
and associated build-up behaviours Pi(t) of the heating 
power are determined, wherein i = l, 2, 3, is an 

index. In a second step a set of calibration curves 8 
of the threshold value time t s is determined as a 
function of the flow rate v for the build-up behaviours 
Pi(t). In a third step, on account of a previously 
measured flow rate or based on a-priori information 
about the presumed flow rate, a preferred calibration 
curve 8 is selected according to a desired measuring 
precision for the flow rate v and according to a 
desired measuring duration t s , and is used to determine 
the flow rate v, or, starting from the calibration 
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curve 8 associated with the lowest flow rate value v i=1 
and rising successively to higher flow rates values v i>1 
or by estimating in a single step, a preferred 
calibration curve 8 is determined according to a 
desired measuring precision for the flow rate v and 
according to a desired measuring duration t S/ and is 
used to determine the flow rate v. 

This gradual procedure is shown in Figure 4 for three 
calibration curves 8, which were obtained with three 
heating power curves Pi(t) at v x =0.25 m/s, P 2 (t) at 
v 2 =l m/s and P 3 (t) at v 3 =2 m/s. A linear measuring 
characteristics t s (v) is given at least pointwise by 
the calibration curves 8. At the same time typical 
values of R s ^35 W/K and Rj*5 W/K were assumed. Even 
with such a disadvantageous size distribution of the 
thermal transfer resistances, namely R S »R I/ the method 
according to invention can still be used to carry out a 
precise and fast measurement in a large flow rate 
measuring range. In principle, very low flow rates 
v<l m/s can still be measured also with the calibration 
curve at v 3 =2 m/s. With the method according to the 
invention, however, the measuring duration t s can ]f 
substantially be shortened and the heat energy : 
requirement can be lowered accordingly. 

Advantageously a number and distribution of the 
calibration curves 8 are selected according to a 
desired measuring resolution and a desired measuring 
range of the flow rate v. 

Figure 1 shows the device for implementing the above- 
mentioned method. Here 2a designates a heating control 
for generating heating pulses 7 for the heating means 
la and 2b designates a measuring device for evaluating 
the thermal measurement and for determining a flow rate 
or flow velocity v or a mass flow from a flow-dependent 
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threshold value time t s until a preset temperature 
threshold value T s is reached at or on the sensor means 
lb. According to the invention the heating control 2b 
comprises means for generating a non-constant heating 
power P with a substantially sublinear build-up 
dynamics P(t) as a function of the time t, and the 
control and evaluating processor unit 2 comprises means 
for at least partial compensation of a nonlinear 
behaviour of the threshold value time t s as a function 
of the flow rate v. 

In a preferred embodiment the control and evaluating 
processor unit 2 comprises hardware and/or software for 
generating a build-up dynmics P(t) proportional to t m 
and/or proportional to a time- independent amplitude 
factor (1+Rs/Rj)" 1 . In addition, the control and 
evaluating processor unit 2 can have calibration means 
2c for performing the first and second procedural steps 
of the above-mentioned calibration procedure. 
Preferably the sensor element 1 comprises an electric 
heating wire la, lb with a temperature -dependent 
resistance, which can be run simultaneously as heating 
means la and as sensor means lb. 

The inventive method and the device for implementing 
the method are suitable for arbitrary fluids 3, in 
particular for liquids 3 or gases 3. 
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Legend 

1 sensor element 
la heating means 
lb sensor means 

2 control and evaluating unit 
2a heating control 

2b measuring device 

2c calibration means 

1/ 2 flow rate sensor, mass flow sensor 

3 fluid; liquid, gas 

4 flow profile 

5 flow channel, pipe 

6 rectangular constant heating pulse (prior 

art) 

7 sublinear non-constant heating pulse 

8 calibration curves 

10 surface of the sensor element 

P heating power, heat output 
P(t), Pi(t) build-up behaviour, increase 

T temperature 

T s threshold value temperature 

T F undisturbed fluid temperature 

t time, time variable 

to/ ti, t 2 times 

t s threshold value time 

v, vi flow rate 



